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Abstract 

The  solution  of  the  initial  value  problem  for  Bogoliubov's  functional 
differential  equation  of  non-equilibri-um  statistical  mechanics  is  obtained. 
This  solution  is  then  expanded  in  an  infinite  power  series  in  the  density 
which  has  the  advantage  that  the  calculation  of  the  leading  terms  requires 
the  solution  of  s-body  problems  only  for  small  values  of  s .   A  derivation 
of  the  equilibrium  equation  by  reduction  from  the  non-equilibri\im  equation 
is  included. 
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1.   Introduction. 

The  statistical  mechanical  treatment  of  a  classical  many-body  system 
usually  begins  with  an  'n-particle  density  function',  D  ,  which  is  the 
solution  of  an  initial  value  problem  for  Liouville's  equation.   There  are, 
however,  two  major  difficulties  with  this  approach: 

1.  In  the  problems  of  interest  the  solution  of  Liouville's  equation 
is  equivalent  to  the  solution  of  an  n-body  problem  where  n  is  very  large, 
and  is  therefore  not  practical. 

2.  The  initial  conditions  are,  in  general,  unknown. 

In  an  attempt  to  circumvent  these  difficulties,  one  introduces  's-particle 
density  functions',  F  ,  defined  by  appropriate  integrals  of  D  .   N.N.  Bogoliubov 
has  shown  '-  -^    that  for  these  functions,  the  Liouville  equation  can  be  replaced 
by  a  fiinctional  differential  equation  for  a  generating  functional  L  [u]  which 
generates  the  functions  F  ,  and  has  obtained  an  expansion  of  the  solution  of 
the  equation  to  first  order  in  the  density. 

In  Section  2  of  this  paper  we  derive  the  functional  differential  equation 
by  a  slight  variation  of  Bogoliubov 's  method.  The  resulting  equation  (20)  differs 
slightly  from,  but  is  equivalent  to,  the  equation  of  Bogoliubov;  however,  the 
form  of  (20)  facilitates  a  new  method  of  solution. 

Section  5  contains  the  main  result  of  this  paper.   In  that  section  we 
obtain  the  solution  of  the  initial  value  problem  for  (20)  by  a  method  similar 
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to  the  method  devised  by  B.  Z\jmino '-  -^    for  the  equilibrium  case.   The  solution 
is  then  expanded  in  an  infinite  power  series  in  the  density.   In  this  fonn, 
it  has  the  advantage  that  for  small  densities  it  may  be  approximated  by  a 
few  terms  of  the  expansion.   Then  to  obtain  an  explicit  expression  for  F 
where  s  is  small,  only  certain  k-body  problems,  where  k  is  small,  need  to  be 
solved.   Furthermore  only  the  Initial  data  for  certain  functions  F.,  where  j 
is  small,  are  required.   If  these  data  are  known,  our  expansion  circumvents 
both  of  the  difficulties  enumerated  above. 


Sections  ^l-  and  5  are  included  for  the  sake  of  completeness .   In  Section 
k   we  carry  out  a  suggestion  of  Zumino  and  derive  the  fxmctional  differential 
equation  for  the  equilibrium  case  by  reduction  from  the  non- equilibrium 
equation.   In  Section  5  we  solve  the  equilibrium  equation  by  a  slight 
simplification  of  Zumino 's  method. 

2.   Derivation  of  the  Functional  Differential  Equation. 

We  consider  a  classical  mechanical  system  of  n  identical  monatomic 
particles  contained  in  a  finite  volume^  V.   The  dynamical  state  of  the  j-th 
particle  is  described  by  the  6  component  vector  x.  =  (q.,  p.)  = 

fa  ,  a  ,  a' ,  D  ,  P  ,  P.)  where  the  q.  are  the  cartesian  coordinates  of  the 

\^y     Hy     Hy     ^y     i^  y     ^y  h^ 

particle^  and  the  p.  are  the  conjugate  momenta,   x .  is  a  point  in  the 
J  J 

phase-space  0      defined  by  the  restriction  that  q.  is  a  point  in  the 
finite  voltmie  V.   The  Hamlltonian  of  the  system  is  given  by 

(1)     H^  =  i:  h(x. )  +  u^, 

i=l 


(2)     h(x.)  =  T(p.)  +  u^(qj.); 


(5)      U^=   H  (^(|q.  -qj), 

l<i<j<n  ^ 

,  a^2 


(^)     T(p.)  =  ^  =  Z 
a=l 


where  m  denotes  the  mass  of  a  particle,  (})  is  the  inter-particle  potential 
and  u^(q.)  is  the  potential  due  to  the  containing  boundary.   Thus  u^(q) 
is  constant  inside  V  and  rapidly  approaches  infinity  at  the  boundary. 
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The  statistical -mechanical  tehavior  of  the  system  is  described  by 
the  n-particle  'probability  density'  function,  D  (t,  x, , . . . ,  x  )  which 
is  symmetric  in  the  variables  (x^,...,  x  ),  is  normalized  by  the  conditi 

(5) 


and   is   a   solution  of  Liouville's   equation^ 


( f:\  n        l-i_L.T^Iv~'\  I  n  n  nn 


Let  s|   denote  the  solution  operator  of  the  n-partlcle  mechanical 
system,  i.e.,  if  the  system  at  time  t  =  0  is  represented  by  the  state 
(x  ,  ...,  X  ]  ,  at  time  t  it  will  be  represented  by  the  state   [x ',...,  x') 

=  si    fx^,...,  X  >  .   Under  suitable  conditions,  the  solution  operator  exists, 
but  of  course  cannot  be  calculated  explicitly  except  when  n  is  very  small. 
If  g  is  a  function  of  ( t,  x  ,...,   x   )  it  is  convenient  to  define  S^  'g 
by  the  equation 


(7)     S^")  g(.,  x^,...,  x^^^)  =  g(.,  s[-)(x^,...,  x^)  ,  x^^^,...,  x^^^), 


In  terms  of  the  solution  operator,  one  may  express  the  solution  of  the 
initial  value  problem  for  Liouville's  equation  in  the  form 


Thus  S^        acts  on  the  first  n  of  the  variables  x.  appearing  in  g. 
t  J 


-  k 


(8)        D  (t,  X,,...,  X  )  =  S^!"^  D  (0,  X,,...,  X  ) 


However,  since  S    cannot  "be  calculated;,  and  since  D  (O,  x  ,  ...,  x  ) 
is  in  general  unknown,  the  solution  (8)  is  of  no  practical  value. 
We  introduce  the  s-particle  density  functions 


(9)       F   (t,  X,,...,  X  )  =  v^   Id 
^^'  n,s^  '   l'    '   s'       j    " 


(t,  x^ ,  .  .  .,  X  )dx   ^  .  .  .dx  ; 
^  •"   l-"    '      n'      s+1     n' 


s  =  0,  1,  2, 


It  follows  that  F    is  symmetric  in  (x  ,  ...,  x  ),  F    =1,  and 

(10)  I     ^F^,3^r--^s=    j     \^r--^n-^'    B  =  l,  2, 


(11)         L^[t,u]  =  J   Djt,  x^,  ...,  x^)  TF   [I  +  vu(x^)]dx^...dx^ 


which  is  defined  on  the  domain  of  functions  u(x)  for  which  the  above 
integral  converges.   By  functional  differentiation  we  obtain 


See  e.g.  [5] 


6  L  s  ,     /  n 

(12)  ^-? V — ^^ r  =  y  "\  ,     D  (t,x,  ,...,x  )  TT    ri+vxi(x.)1dx  ,  .  .  .dx  ; 

^1'     ^    s'  ^         '        I  i=s+l 


•  V 
s  =  0,  1,  2,  .. .,n  ; 

5^L 


(13)  ^-7 \ T"? T      =  — F    t,x,,...,x   ;  s  =  0,  \,    2,...,n. 

^  ^'      5u(x^  )  .  .  .5u(x  )         s,    ,  ,    n,s^  '    V         '    s"  >    +>      >         > 

^1      ^  s'  ]-u=o   n  (n-s)  !     ' 

With  the  aid  of  (13),'  L  may  now  be  expressed  as  a  (finite)  series 
expansion  aroimd  u  =  0: 

(11+)   Lrt,u1=l+y"  ^  (1-  i)...(l-  ^^)  (  F   (t,x,,  ...,x  )u(xj...u(x  )dx,  . 
s=l  '  ' 

A   differential  equation  for  L  may  be  obtained  by  multiplying  (6) 

n 
by  ~[T   lj-+  ^^(x-C  and  integrating  with  respect  to  x  ,  ...,x  over  o 
i=l         ^  1      n       V 

We  obtain 

(15)  ^  =  i:   I   [l+vu(xj^)]  |h(x^);D^  IT  [l+vu(x. )]  I  dx^.  .  .dx^ 
""v  ^  iA        ' 

+  r  f    [l+vu(x    )]  [l+^.Ti(x    )]    J(t)(|q  -q    |);D     fj"        [l+vu(x  Jj  I  d^    .  .  .dx 

l<r<s<n  /  li=l  ijxn 

'V  i^r,s  ' 

By  making  use  of  the  symmetry  of  D   (15)  becomes 


-  =  n  j  [l+vu(x^)]  Jh(x^);  f  D^  f[  [l+vu(x.  )]dX2.  . -dx^  dx^ 
^n(n^|  [i^,^(^^)]  Qi^^^(^^)]  L|q^-q^|;;  |  D^  ft  [:i+.aa(x .  )]dx^ .  .  .dx  l  dx^dx. 

-^2    \    ClWx,):  [l.vu(x,)]  kk.-q^l);  5u(!jgu(x^)  I  ^1^^  • 


We  now  let  n  — >  oo    and  V  — >  oo    in   (l^)    and    (l6)    in   such  a  way  that 


V  =  —  is  finite.      If  we   set 


(17)      Lpt^uJ    =  lim       L    rt,u] 
n— >oo 
V->oo 


then  from  (l4) 


00  /• 

(18)      L[t,u]=l+II     JT     j   F^(t,x^,  ...,x^)u(x^)...u(x^)dx^. 


.dx 


(19)      F^(t,x^,...,x^)    =lim       F^^^(t.x^....x^) 


n->oo 
V->ao 


From   (l6) 


6L 


^^Pi^'5i:u7) 


.||    [;(x^).g[;(x^).g 

It  follows   from  (l8)   that 

(21)   F  (t,x,,...,x  )  =  ^  /   T* ^-7 — X 

^   '   s^  '  l'    '  s'^   6u(x^  ) .  .  .6u(x  ) 


!V^2l^^5u(xj6u(x2) 


Equation  (20)  is  the  functional  differential  equation  which  we  shall 
solve  in  the  next  section.   The  solution  L  is  called  the  'generating  functional' 
because  it  generates  the  functions  F  by  means  of  (21).   We  have  derived  (20) 
by  the  method  of  Bogoliubov  '-  -"  with  a  slight  modification,  and  our  equation 
apparently  differs  slightly  from  the  corresponding  equation  (7.9)  of  [l]  . 
However,,  the  difference  is  only  apparent.   The  two  equations  can  be  shown  to  be 
equivalent,  and  it  will  be  seen  that  our  form  is  more  suggestive  of  how  to 
proceed  in  solving  the  equation. 


By  applying  the  operator 


5u(x^). ..6u(x  ) 
the  infinite  system  of  equations 


dF 


^22)  ^  =  [h,;fJ.  ^ 


l<i<s    ^  ^^-^    ^^-^ 


to  (20)  one  can  obtain 


^.^,'    s  =  1,  2,..., 


(23)   H  =Z  T(p  )  +U  ;  U  =1:     (^(|q-q  I);  s  =  1,  2,  3,...   .    ' 
1=1     ^      ^    "^    l<i<j<s     "-   J 

The  system  (22)  is  equivalent  to  the  single  equation  (20). 

3.   Solution  of  the  Functional  Differential  Equation. 

In  order  to  solve  (20)  we  begin  by  examining  the  case  of  zero  density, 
—  =  0.   We  shall  use  superscripts  'o'  to  denote  this  case.   Thus 

^      1      i  „o 


s=l 


(2J+)   L°[t,w]  =1+^  ij-   F°(t,x^,...,x^)w(x^)...w(xjcbc^...d^^ 


(25) 


^(Pl)^  5^  ^1-  -2    J  Wx,)w(x^)|^(|q^-q^|);^^(j£^^^ 


(26)   F°(t,x,,...,x  )  =  ^-^— AV-7— Y 
^       '        s^  '  1'    '  s'    5w(x  ) . . .8w(x-  ) 


and  (22)  reduces  to 


s  =  1,  2,... 


The  solution  of  (27)  is  immediately  obtained  in  terms  of  the  solutic 
:(s) 


(28)   F°(t,x^,...,x^)  =  sf^)F°(0,x^,...,x^);  s  =  1,  2, . . .    . 

See  e.g.  [ij  . 

It  is  convenient  now  to  denote  the  arbitrary  testing  functions  by  w  instead  of  u. 


Inserting  this  expression  in  (24)  we  obtain  the  solution  of  (25)  subject 
to  the  initial  conditions 

00      r 

(29)  L°[0,w]  =  1+  Z  —   F°(0,x^,...,x  )w(x  )...w(x  )dx  ...djc  , 

s=l   ■  / 

where  the  F  (0,x,  ,...,x  )  are  the  given  initial  data. 

In  order  to  solve  the  general  equation  (20)  we  observe  that  the  form 
of  the  latter  suggests  that  we  try  a  solution  of  the  form 

(30)  L[t,u]  =  L°[t,w];  w(x)  =  u(x)  +  ^  . 

^^^^  M^  =  5^7(^  '  Su(x^)5u(x2)  =  5w(x^)5w(x2)  '   ^  =  W  '      ^^'^  inserting 
in  (20)  we  see  at  once  that  that  equation  is  satisfied  by  virtue  of  the  fact 


But  (20)  must  be  solved  subject  to  the  initial  conditions 

00  / 

(31)  l[0,u]    =X[u]    =  1+  X;;     ^  J   F^(0,x^,  ...,x^)u(x^)...u(x^)d^^...dx^. 

The  functions  F  (0,x  , ...,x  )  are  the  given  initial  data.   In  terms' of  (30) 
this  becomes 

(32)  L°[0,w]  =5r[u];  u  =  W  -i  . 

The  main  result  of  this  paper  is  the  general  s^olutlon  of  (20)  defined  by 
(30).   If  L°[t,w]  is  the  solution  of  the  initial  value  problem  for  (25)  with 
initial  conditions  (52),  then  L[t,u3  is  the  solution  of  the  initial  value 
problem  for  the  general  equation  (20)  with  initial  conditions  (3l)- 

The  method  we  have  used  in  obtaining  this  solution  closely  resembles 

M 

the  method  devised  by  B.  Zumino  •-  -^    to  solve  the  corresponding  functional 
differential  equation  for  the  equilibrium  case.   This  is  discussed  in 
Section  5. 


We  now  proceed  to  obtain  expansions  of  the  fimctions  F  (t,x, ^...,x  ) 
in  powers  of  the  density,  — .   For  this  purpose  it  is  convenient  to  introduce 
a  fiinctional  of  two  variables 


(35)  Q[t,a,w]  =  1*£  iy  I  '^-Vou(J.Lu^)   w(x,)..."(x,)dx,... 


dx^. 


Now  from  {2k)   and  (32) 


6V[0,w] 


\(0,x^,...,x^)   =  5w(x^)...Sw(xj^) 


5u(x^)...5u(x^) 


Hence  from  (28) 
Oh)      Q[t,u,w] 

From  (21),  (30)  and  (3^) 


1=1+1:  |t   F°(t,x^,...Xj^)w(x^)...w(xj^)dx^...dXj^  =  L°[t,^ 
-     k=l    ' 


6^L°[t,w] 
(35)  F^it,x^,...,x^)   =  5^(   )...5^(^  ) 
^1'      ^  s' 

Now,  from  (33) 


^J,        6w(x^)...6w(x^) 


1       1 
w=  — ,  u=  -  — 

v"       V 


r-fi            ^^^ 

oo 

^^^^      6w(x^)...5w(x^) 

-  Z_ 

j  = 

and  from   (31) 

:  ^  ...dx 
s+1    s+j- 


(57)  5u(x^)...6u(x   )  =£^  J  ^nW°'V---"n+J-.s) 


Xu(x.   ^)...u(x.    )dx.   ,  . .  .dx 


This  functional  is  needed  in  the  analysis  in  order  to  avoid  expressions 
involving  divergent  integrals . 


We  now  insert  (37)  in  (56).   The  resulting  double  series  can  be 
rearranged  and  evaluated  for  u  =  -w.   We  obtain 


l!2 


6w(x^) . . .5w(x^) 


k=o 


V:  J.'(k-j).'   -t     k+s^  '  1'    '   k+s' 


^(^s+l)---"^^s+k^^s+l---'^s+k- 


Thus  from  (35) 

oo   1  V  / 
(38)  F^(t,x^,...,x^)  =11  i^f 
k=o     '' 


■^  j'(k-j):   -t    k+s^  •"  l-*    '    k+s' 


s  =  1,  2,  ... 

This  is  our  expansion  of  F  as  a  power  series  in  the  density. 

In  order  to  check  the  results,  one  can  show  in  a  straightforward  manner 
that  (38)  satisfies  the  system  of  equations  (22).  To  verify  that  the  initial 
conditions  are  satisfied,  we  may  set  t  =  0  in  (38)-   Since  S    is  the  identity 
operator,  the  integrand  in  (38)  reduces  to 

(39)  F.   (0,x,,...,x,   )  y~  \,f{    ..,      . 
^■^■^'       k+s^  '    1'        '    k+s'  ^       j'(k-j).' 


But 


(*o)  E  i^^  -iyi:  (J)(-i)'^-J  -  iyd-i)'' =  0,  for.»i,  2,..., 


k   ,  ,  vk-.i     .   k 

z 

thus  every  term  in  (38)  vanishes  except  the  first  and  the  series  reduces  to 
F  (0,x^,...,x  )  as  required. 

The  series  expansion  (38)  is  a  very  useful  form  of  the  solution.   We 


by  terminating  the  series  after  a  few  terms .   Now  the  functions  F  of  main 
interest  are  those  for  which  s  is  small,  and  for  these  functions,  the 
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calculation  of  the  leading  terms  of  the  expansion  requires  a  knowledge  only 

■  wiiexe    Ji    ±ti    tsma_LX    ana.    nnm^ia  i     aax.a     t   ..  ^,^^. ^..  . 

j'     '     1'  '     J' 

where  j  is  small . 


of  solution  operators  S)       where  k  is  small  and  initial  data  F .{0,x^ , . . . ,x .) 


The  leading  terms  of  (38)  are  given  by 


(41)   F^(t,x^,...,x^)=s(^>F^(0,x^,...,x^) 


+  i  /  [:s(-l)F^^^(0,x^,...,x^^^)-s(^'F^^^(0,x^,...,x^^^)]dx^^^.  0(^5) 

'  V 

In  [ij^  Bogoliubov  obtains  the  equation 


(42)   F^(t,x^,...,x^)=s(=)F^(0,x^,...,x^) 


HK-ll 


■^,^,    «"  o(t;) 


1<1<S 


With  a  little  manipulation  it  is  possible  to  reduce  ( k2)   to  the  simpler  form  (^l). 

k .      Derivation  of  the  Equilibrium  Equation. 

In  this  section  we  shall  derive  the  well-known  functional  differential 
equation  for  the  equilibrium  case  by  reduction  from  the  general  equation  (20). 
The  first  step  is  to  derive  a  new  form  of  (20)  by  expanding  the  Poisson  brackets 
that  appear  in  that  equation  (as  is  done  in  (6))  and  by  using  the  following 
identity  which  is  obtained  by  interchanging  integration  variables: 

(45)  I  [u(x,,.  i]  [u(x,).  i]  -^^   -  ,,;^|-^  ax^ex. 


)  b(-i)-  vj  P'^^'^  vj— :;^ —  ^  au(x^)5u(xp  ^1^,- 


With  the  aid  of  (h^) ,    (20)  now  becomes 
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-)i  =  i:l-ME<v^g^ii|i^-i 


n^'Klqi-.ijI)    a         5^ 


Sq^  ^Pj 


Let  us  now  consider  time-independent  solutions, 

CO       , 
(^5)  l[u]  =  1+  Z;  J7   FJx^,  ...,x^)u(x^)...u(x^)dx^...dx^. 
s=l   ■  ' 

Then  ^  =  0,  and  {kk)   will  be  satisfied  if 


'kG)     Yl 


1  a  S     5L 


il  ^^Hqi-^gl)    a         6^1 


1  -^  a  6u(x, 


TI-j  b(-2)-vJ— ^ ^  6u(x^)5u(x^)  ^2 


dx^^  =  0. 


M{  ^V, 


This  equation  is  an  identity  in  x  =  (q  ,p  ).   It  is  sufficient  for  {kk) , 
but  not  necessary. 

Following  a  suggestion  of  Zumino '- ^ ,  let  us  now  consider  solutions  of 
(46)  for  which 


{kl)      F^(x^, ...,x^)  =  c"   exp 


where  6  is  a  constant  and 


1    r2  2"! 


^s^V-'-^s) 


(i+8)      c   =        exp 


L'   2meJ 


dp. 

Let  l[uJ   denote  the   restriction  of  L  [u]   to  the  domain  of  functions  u  =  u(q) 
which  are   independent   of  p.      Then  if   (^7)    is  assumed. 


u(x^)cbc^...dx^ 


00       -,        /  k  -p 

and 

(50)      L[u]    =1+YL     tr        fj,(qi.---.qi,)u(qi)...u(q3^)dq^...dq 


13 


By  functional  differentiation  of  (kg)   and  (50)  it  is  easy  to  show  that 


(51) 


6"L 


5u(x^)...5u(Xg) 


=  c    M   exp 
i=l 


2rae 


5°L   

6u(q^)...&u(qg) 


and 


(52) 


a  6u(x  )5u(x  )     me  5u(x  )6u(x2) 


With  the  aid  of  (52),  {^6)   hecomes 


6u(x^)Bu(x2)  "^2 


If  we  restrict  (55)  to  functions  u  =  u(q)  and  use  (51)  we  obtain 
2  n,     f 


(5^)   c'^exp 


-1 
c  exp 


2me 


dp. 


dqoV  =  0, 


(55) 


a    SL    .  1 


[-(%Ha 


.S<^(IV^2l)      6^ 


^BR^-^eJL-^^2^-vJ— ^^^^       6u(q^)5u(q2)  ^^2 


dq^  =  0;  a  =  1,  2,  5 


(55)  is  the  well-known  equation  of  equilihrliun  theory,  where  0  =  kT,  k  is 
the  Boltzman  constant,  and  T  is  the  absolute  temperature.  This  equation  is  usually 
derived  from  an  assumption  about  the  explicit  form  of  D  .   This  form  is  given  by 


See  e.g.  [l]  or  [k]  . 


-  ih  - 

(56)  D^  =   Z;^  exp[-  I  hJ    •>    \  =  t^  T(p.)+  U^   ;    where 

(57)  Z^=(      exp[-|Hj]dx^..    d.^=0^c-;    Q^=   jexp[-|ujdq^...dc^. 

%  V 

The  purpose  of  this  section  has  been  to  show  that  the  equilibrium  equation 
(55)  can  be  derived  from  the  general  equation  (20)  by  using  the  assumption  (^7). 
It  is  not  surprising  that  this  can  be  done,  in  view  of  the  fact  that  (^7)  is  a 
consequence  of  (56).   To  see  this,  we  use  (19)  and  (9)  to  obtain 


(58)   F^(t,x^,...,x^)  =lim  F    =lim  V^j   Z^^  exp  [-  |  H^]dx^^^.  .  .dx^ 
n->co  '  n->oo    •'         ^  -^ 

V->oo       V->oo    V 

=  lim  V^c-exp[-  ^(P^...■^P^^  j  Q;'  exp[-  |  U^]dq^ 


n— >oo 
V->OD 


From  (58)  we  see  at  once  that  (^7)  follows  with 

(59)  fjq,,...,q3)  =  lim    V^  I  Q;lexp[-  i  U  Idq^^^  .  .  .dq^. 

n->oo    '  •-     -' 

V->oo 

Before  proceeding  to  the  solution  of  (55)  we  point  out  that  that  equation 
is  also  equivalent  to  an  infinite  system  of  equations,  given  by 

^^k   1  ^^k      1   /  ^^^l'll"'lk+ll^ 

(60)  -^  +  1-^  f^-^h 
^  A  a   G  5,  a  k   Gv 

oq-L     dq^ 


See  [1]  for  the  derivation. 


15 


5.      Solution  of  the   Equllibriiun  Equation. 

In  this   section  we  shall   solve  the  equilibrium  equation  (55)  by  a   slight 
simplification  of  a  method  due  to  B.    Zumino  ^  -I.      As   in  Section  5  we  begin  by 
examining  the   case  of   zero  density,   -  =  0.      We  again  use   superscripts    'o'   to 
denote  this   case  .      Thus 


(61)      L°[w]=l+£     ^    lfl{^^,...,q^)wi<i^)...wU^)dq^...dci^ 


(63)  flU^,.--,%)    =  5w(q^)...6w(qJ 
and  (60)  reduces  to 

(fi"*'  Zi^-.zi'l"'-'      «  =  1,  2,  3,  k  =  l,  2,...    . 

In  order  to  solve  (Sk),    set 

(65)  f°=  Cj^(q^,...,q^)exp[-|uJ  ;  k  =  l,  2,...   . 

From  (64), 

(66)  Jk  .  0  ;  a  =  1,  2,  3;  k  =  l,  2,...   . 

Since  C,  (q^ ,  . .  .,  q,  )  is  symmetric  in  its  argiaments,  it  follows  that  C,  is  a 
constant.   In  order  to  determine  the  constant,  we  observe  that  by  letting 
n->oo,  V->a3  in  (lo)  we  obtain 
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(67)   llm   ^   M's  "^l- 


V->ao    V 

Now  from  (^7) 

(68)     lim       ^ 
V->oo   V 

V 

and  from   (65) 

(69)      lim       — 
V->oo   V^ 

exp 

V 

dq^  .  .  .dq^ 


[-|uj  dq^...dq^.i 


It  is  clear  from  (69)  that  C  =  1  for  potentials  (j)(r)  which  vanish 
sufficiently  rapidly  as  i^>oo .   We  shall  therefore  Impose  as  a  condition  on 
(j)  that  the  left  side  of  (69)  he  equal  to  1  for  s  =  1,  2,  5;,  •  •  •  •   It  follows 
now  from  (65)  that 

(70)  f°(q^,...,q^)  =  exp[-|u^]  ;  k  =  l,  2,...   . 

The  solution  of  (55)  i'or  non-zero  density  can  be  obtained  from  the  zero 
density  solution  in  a  manner  very  similar  to  the  procedure  used  in  the  non- 
equilibrium  case.  We  begin  with  a  trial  form  of  the  solution  slightly  more 
general  than  the  one  used  in  Section  3: 

(71)  L[u]  =  L°[w];   w  =  a(u  +  -) ;    a  =  const. 
By  functional  differentiation  we  have 


(72) 


5u(q^)...5u(q^)     ^   6w(q^) . . .6w(q^) 
and  substituting  in  (55)  we  see  that  the  latter  equation  is  satisfied  because 
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In  order  to  determine  the  constant,  a,  we  observe  first  that  since  U 
is  a  function  only  of  the  coordinate  differeftces  (q. -q.),  (59)  implies  that 
f  (q  )  is  a  constant,  and  (68)  implies  that  the  constant  is  unity.   Thus 


(73)  fi(qi)  =  1. 

Now  from  (73)  and  (50)  it  follows  that 


(7^)  ^-Fh     =  1;  '^[0]  =  1- 

^'     '      5u(  q  )   I  u=o    '    L  -J 


This  in  turn  implies  that 

We  shall  see  that  (75)  siof flees  to  determine  the  constant  a. 
Now  from  (50)  and  (72), 


6°L 


S-^L^ 


(76)   f^(q^,...,q3)  =  5^(q^)...5^(^^)   ^=0=^   6w(q^) . . .5w(q^) 

It  would  appear  that  we  need  only  differentiate  (61)  s  times  and  set  w  =  — 
to  obtain  an  explicit  formula  for  f  (q  , ...,q  ).   However,  this  is  incorrect 
because  f  ^»  1  for  large  |q.-q.|  and  the  Integrals  in  (61)  converge  only  for 
testing  functions  w(q)  which  vanish  sufficiently  rapidly  at  infinity.   For 
w  =  — ,  the  Integrals  diverge.  The  difficulty  is  that  (61)  does  not  represent 
the  functional  L  [w]  in  a  sufficiently  large  domain  of  fiinctions  w(q).   What  is 
needed  is  an  'analytic  continuation'  of  the  representation  of  the  functional. 

Such  an  analytic  continuation  can  be  obtained  by  the  following  transformation 

M 

which  was  suggested  by  Zumino  '-  -'  . 
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(77)      L^    H    =   exp  M°[wJ      ,    M°[w]    =  log  L°  [w]    , 


_-  CO  I 

(78)      M^[w]=^     ^        g°(q^,...,qj^)w(q^)...w(qj^)dq^...dqj^. 


„0/       V  5L 


|m°M 


5M° 


&w(q-^) 


5w(q^ ) 
w=o  ^  ^1 ' 


-.^%). 


Proceeding  in  this  manner,  we  may  obtain  the  following  relations  between  the 

(79)  f°(qi)  =  g°(qi) 

f°(q^,q2)  =  g°(q-[_)g°(q.2)  +  S2(li'^2^ 

^°^i\>^2''^^^   =  §1(^1)6^(^2)61(13)  +  g2(<li^q2)®?S)  "^  g2^^2''^3)%^^l) 
+  g2(li^q3)g°('l2)  "^  65(11^  ^2'%)' 

(80)  g°(q^)  =  f°(q^) 

g2(v^)  =  ^l^^r^2^    -   f  1(11)^1(^2) 

g3(qi,  V^j)  =  f3(V  V^j)  ■  ^l^^l)^2(  Vlj)  -  ^\^<^)^\^^^,^^)    -   f\W^)fl{^^,^) 
+  2f°(q^)f°(q2)f°(q5)  , 


large  values  of  |q. -q. | 
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Let 

(81)      z  =  f   . 

From   (77)    and  (75), 

fOr.\                     5M 

1 

sfP 

1 

^^^^       5w(q^) 

^=^  ^[^: 

6v(ci^) 

a 
w=z 

and  from  (76),  (77),  and  (75) 


(85)   f2(qi,q2)  =  a  ^^^^— ^^^^  I  ^^  =  a  L  [z]   ^^^  ^^^  +  5w(q^)5,(q^) 


52^?^ 


2  ll_     ,  5^^?^ 

^2  ^   6w(q^)5w(q2) 


Thus  from  (78) 


1   k    o 


(8M   f2(q,,ci2)  =  1  +  a^  i:  t^      ]    W^r ---W^^^J-'-^W' 
k=o       ^ 

and  from  (82) 

(85)  1  +E  ^  z""  I  g°^i(q,v-^k^^'ir--^^k  =  a  • 

k=l       / 


By  virtue  of  the  remark  made  at  the  end  of  the  last  paragraph,  we  see  that  the 
integrals  appearing  in  (84)  and  (85)  are  convergent.   If  ve  set 


(S6)  \  =  h\  4^^'^i'---'\-i^^\---^\-v  ^  =  2.  3, 


..;  b  =  1; 


then  (85)  takes  the  form 

(87)   V  £  kb  z^  =  1. 
k=l   ^ 

The  b  are  called  'cluster  integrals'  and  are  independent  of  q  because  the 

f°,  and  hence  the  g°,  are  f\inctions  only  of  the  coordinate  differences.  The 

quantity  z  is  called  the  'activity'.   (84)  expresses  f^  as  a  pover  series  in  the 
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ictivity,  and  the  latter  is  related  to  the  density  —  by  (87).   In  order  to 
)'btain  an  express 
CO   a. 

(88)  z  =  j:  ^ 

j=l  v'J 
and  insert  in  (87).   One  obtains  easily  that 

1   2b 

(89)  2.  =  ^  -  ^  +  O(^)  , 

V  V 

and  inserting  in  (8k)   we  obtain 

(90)  f^{q^,1^)    =   1+  g2(q^,q2)+  ^[j  g°(  q^,  Og' ^^J  ^"^^^'^^^  ^^l' ^^2^  j  ^^^l' ^J^^S]""  °^~2^ 
It  can  easily  be  shown  that  this  result  agrees  with  previously  given  expressions 


order  — .  Formulas  for  f  for  s  >  2  can  be  obtained  by  an  obvious  generalisation 
However,  the  solution  of  the  functional  differential 
equation,  (55),  is  in  principle  already  given  by  (71)^  where  L  is  given  by  (6I) 
and  (70),  and  a  is  determined  by  (87). 
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